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Abstract 

Using a variation of Lueschers geometric charge definition for SU(2) lattice gauge 
theory, we have managed to give a geometric expression for it's Chern-Simons term. 
From this definition we have checked the periodic structure. We determined the 
Chern-Simons density for lattices L 4 and L 3 x 2, 4 with L = 4, 6, and 8 near the 
critical region in the SU(2) Higgs model. The data indicate that tunneling is in- 
creased at high temperature. 
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1 Introduction 



Quite some time ago 't Hooft found that the baryon number and the lepton number 
are not conserved in the electroweak theory Jl] B — L is of course conserved due 
to the anomaly cancellation. Baryon number violation is caused by the nontrivial 
topological winding of the SU(2) gauge fields. The anomaly of the fermionic current 
relates the winding of the gauge fields and the change in the baryon number by an 
amount 

B(t 2 ) - B{ti) = ^ jf * / (PxtriF^F^] (1.1) 

where Nf is the number of families of quarks and leptons. Equivalent, in the tem- 
poral gauge A = we can relate the change in the baryon number to the change 
in the Chern-Simons number 

B(t 2 ) - B(h) = N f [N CS (t 2 ) - N cs (ti)] (1.2) 

where the Chern-Simons number Ncs is 

N cs = -g^ / ^xe^trlMdjAk + ^AjA k )]. (1.3) 

At zero temperature such processes are exponentially suppressed as exp( — 2n/aw), 
a ~ 1/30. This is because any gauge field configuration which changes the winding 
number has an action at least that of the barrier height 2ir/aw- 

At high temperatures which prevail in the early universe such an exponential 
suppression is absent since the system can pass over the barrier classically. The 
only suppression factor is the Boltzmann factor exp(—/3E) where E is the barrier 
height, and it was shown by the authors in ref. || that this factor goes to one. 
Thus, any baryon asymmetry generated at the GUT scale will get washed out as 
the universe approaches the electroweak phase transition from above. This happens 
typically around 10 Tev and is caused by static objects called sphalerons. These are 
unstable, but finite energy solutions of the classical Yang-Mills Higgs fields. While 
an instanton tunnel between two vacua, a sphaleron moves from one top of the 
barrier to the next. If we assume that the vacuum has Ncs — 0, then the sphaleron 
has a baryonic charge of 1/2. 

The radical solution to the problem of the observed baryon asymmetry in the 
universe was put forward by the same authors using CP violating processes in the 
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electroweak theory and assuming thermal non-equilibrium provided by the expansion 
of the universe. In order for the whole scenario to work, it is necessary that the 
electroweak phase transition is of first order, which severely restricts the Higgs and 
the top quark masses. For more on the origin of the baryon asymmetry see ref . || . 

There has been some lattice studies of baryon number violating processes in the 
two dimensional Abelian Higgs model, as well as in the four dimensional SU(2) Higgs 
model ref. || and ref. H . The configurations are prepared at high temperature and 
the system is allowed either to change via the classical Hamiltonian equation of 
motion or by a Langevin equation. It is then possible to study how the Chern- 
Simons term Ncs changes during the evolution as a function of the temperature. It 
was shown that when ANcs = ±1 the system passes through a sphaleron transition. 

Since baryon number violating processes are of non-perturbative origin it is not 
clear if the sphaleron approximation is sufficient enough. There might be other 
relevant configurations contributing. This warrants a study in the full Higgs model 
without assuming high temperature. The quantity we have studied is the constraint 
effective potential V(Nqs) 

exp(-V(N cs )) = JdAd<f>exp(S{A,<l>)) 

S{N CS + 3-2/ (Pxei jk tr[Ai{djAk + (1.4) 

As a measure for baryon number violating processes we can use the parameter 
r = exp(V(0) — V(^)). To evaluate N C s we have used a geometric definition given 
in ref. ||. After vectorizing the code first made by one of us (MLL) , see ref. @, we 
are able to calculate Ncs efficiently. 

2 Topological charge and the Chern-Simons term 
in the continuum 

As a warm up, let us first consider the two dimensional Abelian gauge theory. 
We always work in Euclidean space time. The theory is determined from the gauge 
potentials A^(x), (// = 0, 1) which under a local gauge transformation g(x) transform 
as: 

A^x) -> A^x) + ig-\x)d»g(x). (2.1) 
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The gauge field tensor F^x): 

F^(x) = d,A v {x) - d„A„(x) (2.2) 

is of course gauge invariant. From this we can define a gauge invariant and integer 
topological charge Q: 

Q = ±- [ d 2 xe^F^(x) e Z (2.3) 

with M being the manifold. The topological charge density q(x) can be written as 
a perfect derivative 

q(x) = —e^F^x) = d^K^x) (2.4) 
where we have introduced the Chern-Simons density K^{x). 

= A v {x). (2.5) 

It transforms like: 

K^x) -> K„(x) + ^ v g- 1 {x)d v g{x). (2.6) 
Now we can define the gauge dependent Chern-Simons number Ncs 



N cs = / d Xl K (x) £ Z (2.7) 

JdM=S 1 

Here, the boundary of the manifold dM is assumed to be a one sphere S 1 . In general 
Ncs is n ot an integer. This happens only if the the gauge field is pure gauge, that 
is A^x) = ig~ 1 (x)d tM g(x). However the gauge variation is always an integer 

SN CS = ±-f dx 1 g- 1 (x)d 1 g(x) G Z. (2.8) 

This follows from homotopy theory by considering the mapping 

9(x) : S l - U{1) = S\ (2.9) 

Such mappings are characterized with the homotopy class 

TiiiS 1 ) e z. (2.10) 

Now we are prepared for the four dimensional SU(2) gauge theory. The gauge 
fields are now non-abelian A fl (x) = A^(:r)T a (// = 0, 1,2,3 and T a ,a = 1,2,3 are 
the Pauli-matrices) , and they transform under a local gauge transformation g(x) as: 

A^x) - g-\x)A,{x)g{x) + g-\x)d lx g{x). (2.11) 
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The gauge field tensor F pv (x): 

F^(x) = d^A^x) - d„A^(x) + [A^x), A v {x)\ (2.12) 

transforms now gauge covariantly 

F^(x) -> g-\x)F flu (x)g(x) (2.13) 

In analogy with two dimensions we have a gauge invariant and integer topological 
charge Q: 

Q = / d 4 xe pi/pa tr[F pu (x)F pa (x)] G Z (2.14) 

OZ7T JM 

The topological charge density q(x) is again a perfect derivative 

= -^r e ^P<rtr[F^(x)F pa (x)} = d^K^x) (2.15) 

OZ7T 

where the Chern-Simons density K p (x) is 

K n( x ) = -^^u P atr[A u (x)(d p A a (x) + -A p {x)A a {x))\. (2.16) 
It transforms like: 

K p (x) -> - ^^e^p (T tr^~ 1 (a;)^(a;)^- 1 (a;)9 p ^(a;)^" 1 (a;)9 (T ^(x)] 

'" e AH/ ^r[9 I ,(9^(a;)y- 1 (x)A T (x))]. (2.17) 
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Finally the Chern-Simons number N C s is: 

iV C 5= / d 3 xK (x)(£Z (2.18) 

This time the boundary of the manifold DM is assumed to be a three sphere S* 3 . 
While Ncs is only an integer for pure gauge configurations, the gauge variation is 
an integer (the boundary term vanishes) 

SN CS = -777"^ / d 3 xe 0l/pa tr[g~ 1 (x)d l/ g(x)g~ 1 (x)d p g(x)g~ 1 (x)d a g(x)} G Z 

Z47T JdM=S :i 

(2.19) 

This follows also from homotopy theory using the mapping 

g(x) : S 3 -> SU{2) = S 3 . (2.20) 
Such mappings are characterized with the homotopy class 

n 3 (s 3 ) G Z. (2.21) 



3 Topological charge and the Chern-Simons term 
on the lattice 



We will now consider the lattice versions of the topological charge and the Chern- 
Simons number. We begin with the two dimensional U(l) theory and we assume 
that the manifold is a two torus M = T 2 and we shall cover M with cells (plaquettes) 
c(n). Let the gauge potential A!l~^ be defined on c(n — fi) and likewise A™ be defined 
on c(n). At the faces (links) f(n,fi) = c(n — p) n c(n), we relate the two potentials 
via a gauge transformation or transition function v n ^ 

A n u -* = A n u +tv^d u v n ,,. (3.1) 

Starting from eqn. 2.3 we find for the topological charge 

Q = ^- Yl e ^ L, . dx v v-^ d v v ntli e Z. (3.2) 

Z7r n,n • / /(™>W 

In the local axial gauge in c(n) one has at the corners 

v n ^(x) =w n _ [l (x)w~ 1 (x) (3.3) 

where w n (x) is a parallel transporter. To get to the lattice result one must interpolate 
v n ,ft( x ) continuosly between v n>fi (n) and v nifl (n + v) say 

VnA X ) = KA n + v) v w( n )] Xv nA n ) < X < 1. (3.4) 

After some simple algebra one arrives at: 

Q = 7T S ar 9 [UnfiU^U-^U-^] eZ, -n < arg < tt. (3.5) 

Here U nifl are the links and we notice that the charge only depends on the plaquette 
angle, hence it is gauge invariant. If we write U n ^ = exp[ik n>fi ] then the natural 
definition of Chern-Simons term is 

Ncs = ^Y.Ko (3-6) 

where the sum is over the spatial lattice only. Since under a gauge transformation 
g(n) = exp[i0(n)] 

k n ,o -> k n ,o + 0(n) - 6{n + 6) + 2tt x to, to e Z, (3.7) 
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it is clear that Ncs only can change by an integer. 

This generalizes to SU(2) in four dimensions M = T 4 . The cells are now hy- 
percubes and the faces f(n,fi) are cubes. There are also plaquettes p(n,fi,u) = 
f(n)nf(n — u). For Q one has, see ref. |J. 
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fj,upa 



^ r tVM( a: )^ t '™.M( a; ) u ni(^)^n, M (a;)w„i(a;)5 ff t; niM (x)] 

+3 / d 2 xtr[v n ^{x)d p v-^{x)v~ l _^ v (x)d a v n ^^ u (x)]}. (3.8) 

At this point we shall deviate from Lueschers charge and use a slightly different 
procedure given by Seiberg. It has the advantage that one can define a Chern- 
Simons term also. The topological charge is: 

Q = E?( n ) 

n 

?( n ) = ^EMn^-W^) -7r<g(n)<7r. (3.9) 
Here is the local Chern-Simons term which leads to 

Ncs = ^~J2 k n,o (3-10) 
with the summation over the spatial lattice only. Explicitly: 



d 3 tr[S(x)d u S(x)- 1 S{x)d p S{x)~ 1 S{x)d a S(x)- 1 ] 

f(n,n) 

+3 / d 2 xtr[p- 1 {x)d p P{x)S- 1 {x)d a S{x)]]. (3.11) 
The function P is defined on a plaquette with corners a, b, c, d: 

P(x, y) = UyjUl(U ac U cd U db U ba )yU ab U* d r, < (x, y) < 1. (3.12) 
In particular at the corners of the plaquette: 



P(0,0) 


= l 


P(1,0) 


= E/a6 


P(0,1) 


= f/ac 




= U ac U c d 



(3.13) 



Therefore, P can be interpreted as the gauge transformation on the plaquette which 
brings the links into the complete axial gauge. Likewise S is defined on the cube, 
such that it brings the links into the complete axial gauge. We have not given the 
expressions for S since they are rather cumbersome. Letting U n ^ = exp[aA ntfl ], and 
taking the naive continuum limit a — > one finds: 

la 3 2 

~^~kn,u(x) — — ~€uupcTtr[A ni ,(dpA n(T + —A np A n(T )\ (3-14) 

Z7T 07T Z O 

which agrees with the continuum expression eqn. 2.16. The only difference between 
Luescher's and Seiberg's definition is that Luescher gauge fixes to the complete axial 
gauge in the whole hypercube. All links in the hypercube are therefore typically of 
the form U ac U c dUdbUb a - Thus the k n>lx is by itself gauge invariant and it contributes 
0(a 4 ) in the naive continuum limit. So it can not be interpreted as a Chern-Simons 
term. Also one can relax the restriction of the topological charge given in eqn. 3.8. 



4 Tests and Monte-Carlo results for the Chern- 
Simons density 

We have used an SU(2) Higgs model with action: 

S = -f E tr[U n ,, w ] - Kj2tr[& n U n ^ n+ »] + A($t$ n - l) 2 . (4.1) 

n,p<u p, 

We always tried to work close to the Higgs phase transition, so the gauge Higgs 
couplings were chosen with: A) = (2.25, 0.5) and k = 0.25, 0.30. 

To do the integrals in eqn. 3.10 we have used a vectorized version of the topology 
program used in ref. j7|. One integration is done in analytic form and we are left 
with a two dimensional integral. We have used the following strategy, which turned 
out to be the most efficient. Perform a Gaussian integration with 8x8 points and 
store the results for the eight k n>ll in each cube. Redo the same thing with 16 x 16 
points and compare the results for each k n fl . If the difference is less than 0.001 
we accept the contribution. Otherwise we collect the cubes who's integrals have 
not yet converged and redo these with 32 x 32 points instead. Compare with the 
previous values. Usually, at this point only a few integrals have not converged. 
Those are normally quite tricky, so for these we use a library integration routine 
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with interval adaption. The typically time for one topological charge on a 6 4 lattice 
is 100 seconds on the CRAY-YMP. The charges are integers up to errors of the 
order 10 -4 . To make the Seiberg charge converge, it is necessary to perform a global 
Landau gauge fixing before the integration. This is allowed since the charge is gauge 
invariant. We have checked that Lueschers and Seibergs charge definitions agree in 
each hypercube up to an integer. Atmost a few of hypercubes have a charge outside 
the interval ] — 1/2, l/2[, so that often the two definitions agree. Now for the Chern- 
Simons term we only need to evaluate k n $ for one timeslice. That takes around 5 
seconds. As we know the Chern-Simons term is gauge dependent, but for a check 
we looked at the periodic structure for a 4 3 x 2 lattice. We did 100 configurations 
and we first measured Ncs without any gauge fixing. This required a much higher 
accuracy for our integrals before convergence. In Fig. 1 we have plotted the Chern- 
Simons density. Notice that there are many configurations which have Ncs close 
to an integer. These configurations can be interpreted as being pure gauge. In 
Fig. 2 we have plotted the same quantity but with one Landau gauge fixing sweep. 
There are now more configurations with Ncs centered around 0, ±1. Most of the 
configurations have changed Chern-Simons number by an integer. In Fig. 3 we have 
performed 5 Landau gauge fixing sweeps and Ncs is now centered at 0. Since we 
have demonstrated that the Chern-Simons term only changes by an integer under 
gauge transformations, we can safely restrict N C s to the interval ]-l/2,l/2[. We 
now compared the density for the two lattices 6 4 and 6 3 x 2. For these and all 
the other lattices we used 50 Landau gauge fixing sweeps, this is sufficient for the 
integrals to converge. Typically we have around 1000 Chern-Simons numbers. See 
Fig. 4 and Fig. 5. There is a trend in the direction of a flatter distribution at finite 
temperature. We would like to interprete this as the system likes to tunnel more 
often. We have also done quite a number of other symmetric lattices 4 4 , 6 4 and 8 4 , as 
well as asymmetric 8 3 x 2,4. We emphasize that these results are only preliminary 
and a lot has to be done before we can really say that tunneling is improved at high 
temperature. We also need to see how the density depends on the spatial volume. 
Finally we would like to mention that an alternative definition of the Chern-Simons 
term has been derived by the authors in ref. ||. 



9 



References 

[1] G. 't Hooft, Phys. Rev. Lett. 37 (1976) 8; Phys. Rev. D14 (1976) 3432. 

[2] V. A. Kuzmin, V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B155 
(1985) 36. 

[3] E. W. Kolb and M. S. Turner, The Early Universe, Frontiers in Physics, 
Addison- Wesley Publ. Co; E. W. Kolb talk presented at the HLRZ workshop 
on Dynamics of First Order Phase Transitions, June 1-3, 1992. 

[4] D. Yu. Grigoriev, V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B216 

(1989) 172; Nucl. Phys. B326 (1989) 737. 

[5] J. Ambj0rn, M. L. Laursen and M. E. Shaposhnikov, Phys. Lett. B179 (1987) 
757; Nucl. Phys. B353 (1989) 483; 

J. Ambj0rn, T. Askgaard, H. Porter and M. E. Shaposhnikov, Phys. Lett. B244 

(1990) 479; 

A. I. Bochkarev and Ph. de Forcrand, Phys. Rev. D44 (1991) 519. 

[6] M. Liischer, Commun. Math. Phys.85 (1982) 39; 
N. Seiberg, Phys. Lett. B148 (1984) 456. 

[7] I. A. Fox, J. P. Gilchrist, M. L. Laursen and G. Schierholz, Phys. Rev. Lett. 54 
(1985) 749. 

[8] M. Gockeler, A. S. Kronfeld, G. Schierholz and , U. -J. Wiese, work in progress. 



10 



5 Figure Caption 



Figure 1. The Chern-Simons density as a function of the Chern-Simons number 
Nqs with zero gauge fixing sweeps. The volume is denoted V, and /3, k, A are the 
gauge Higgs couplings. Figure 3. The Chern-Simons density as a function of the 
Chern-Simons number Nqs with one Landau gauge fixing sweep. Same parameters 
as before. Figure 3. The Chern-Simons density as a function of the Chern-Simons 
number N C s with five Landau gauge fixing sweeps. Same parameters as before. 
Figure 3. The Chern-Simons density at zero temperature as a function of the re- 
stricted Chern-Simons number Ncs with fifty Landau gauge fixing sweeps. Same 
parameters as before. Figure 5. The Chern-Simons density at finite temperature 
as a function of the restricted Chern-Simons number N C s with fifty Landau gauge 
fixing sweeps. Same parameters as before. 
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